Critical velocity for superfluid flow across the BEC-BCS crossover 
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Critical velocities have been observed in an ultracold superfluid Fermi gas throughout the BEC- 
BCS crossover. A pronounced peak of the critical velocity at unitarity demonstrates that superflu- 
idity is most robust for resonant atomic interactions. Critical velocities were determined from the 
abrupt onset of dissipation when the velocity of a moving one dimensional optical lattice was varied. 
The dependence of the critical velocity on lattice depth and on the inhomogeneous density profile 
was studied. 

PACS numbers: 03.75.Kk, 03.75.Lm, 03.75.Ss 
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The recent realization of the BEC-BCS crossover in ul- 
tracold atomic gases [l[ allows one to study how bosonic 
superfluidity transforms into fermionic superfluidity. The 
critical velocity for superfluid flow is determined by the 
low-lying excitations of the superfluid. For weakly bound 
fermions, the (Landau) critical velocity is proportional to 
the binding energy of the pairs, which increases mono- 
tonically along the crossover into the BEC regime. How- 
ever, the speed of sound, which sets the critical velocity 
for phonon excitations, is almost constant in the BCS 
regime, but then decreases monotonically on the BEC 
side, since the strongly bound molecules are weakly inter- 
acting. At the BEC-BCS crossover, one expects a rather 
narrow transition from a region where excitation of sound 
limits superfluid flow, to a region where pair breaking 
dominates. In this transition region, the critical velocity 
is predicted to reach a maximum [1, 0, 0] . This makes the 
critical velocity one of the few quantities which show a 
pronounced peak across the BEC-BCS crossover in con- 
trast to the chemical potential, the transition tempera- 
ture Q , the speed of sound 0, 0] and the frequencies of 
shape oscillations which all vary monotonically. 

In this paper, we report the first study of critical ve- 
locities across the BEC-BCS crossover, where a Feshbach 
resonance allows the magnetic tuning of the atomic in- 
teractions, and find that superfluid flow is most robust 
near the resonance. Our observation of a pronounced 
maximum of the critical velocity is in agreement with 
the predicted crossover between the two different mech- 
anisms for dissipation. 

Critical velocities have been determined before in 
atomic BECs perturbed by a stirring beam 
as well as by a ID moving optical lattice [lH. In both 
cases, the inhomogeneous density of the harmonically 
trapped sample had to be carefully accounted for in or- 
der to make quantitative comparisons to theory. Here 



* Present address: Department of Chemistry and Chemical Biology, 
Harvard University, Cambridge, Massachusetts 02138 
Website: cua. mit . edu /ketterle_group 



we mitigate this problem by probing only the central re- 
gion of our sample with a tightly focused moving lattice 
formed from two intersecting laser beams. For decreasing 
lattice depths, the critical velocity increases and, at very 
small depths, approaches a value which is in agreement 
with theoretical predictions. 

In our experiments, we first create a superfluid of ^Li 
pairs according to the procedure described in previous 
work [13]. Forced evaporative cooling of an even mix- 
ture of the two lowest hyperfine states is performed at 
a magnetic field of 822 G, on the BEC side of a broad 
Feshbach resonance centered at Bq = 834 G. This results 
in a nearly pure Bose-Einstein condensate of N ^ 5x 10^ 
pairs in a cross optical dipole trap with harmonic trap- 
ping frequencies i^x,y,z = (65,45,50) Hz. The Fermi en- 
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ergy of the system is Ep — hv {6N) ' =h x 7.6 kHz. To 
form the moving lattice, we focus two phase- locked 1064 
nm laser beams to intersect at the sample with an angle 
of ~ 90° (see Figure 1). The resulting ID lattice has a 
spatial period of Al = 0.75 /im. A frequency difference 
between the two beams of /S.v causes the lattice to move 
with velocity vl = Ai^. The beams have waists of 
20 /im and 60 jam respectively, and address a relatively 
homogeneous region at the center of the cloud which has 
Thomas-Fermi radii Rx,y,z = (63,91,82) /im. The mini- 
mum density at the position of the waist is 42% of 
the central density. 

The lattice which necessarily varies in depth across the 
sample, is characterized by its peak depth Vb specified 
in units of Ep or the recoil energy Er = K^/(^m\\) = 
h X 7.3 kHz, where m is the molecular mass. The lattice 
depth is calibrated using Kapitza-Dirac scattering. Due 
to the inhomogeneity of the lattice, the uncertainty is 
40%. The lattice depths explored in this work are suffi- 
ciently small such that motion induced in the laboratory 
frame is negligible, in contrast to [l3 |. 

The lattice moving at a constant velocity is adiabat- 
ically ramped up and held for a time t up to 2 s, after 
which the lattice is ramped down and all confinement is 
switched off. As in previous work [l^, a fast magnetic 
field ramp is used to reduce strong interactions in order 
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FIG. 1: Onset of dissipation for superfluid fermions in a mov- 
ing optical lattice, (inset) Schematic of the experiment in 
which two intersecting laser beams produced a moving optical 
lattice at the center of an optically trapped cloud (trapping 
beams not shown). Number of fermion pairs which remained 
in the condensate Nc after being subjected to a Vb = 0.2 Ef 
deep optical lattice for 500 ms, moving with velocity vl, Sit 
a magnetic field of 822 G (l/krCL = 0.15)) An abrupt onset 
of dissipation occurred above a critical velocity Vc, which we 
identify from a fit to Equation 1. 



to probe the center-of-mass momentum distribution of 
the pairs. Subsequently, absorption imaging is done on 
the atomic resonance line at 730 G. A bimodal fit re- 
veals the number of pairs remaining in the condensate 
A/'c, providing a measure of the heating incurred during 
application of the moving lattice. 

Figure [T] illustrates the characteristic dependence of 
dissipation on the velocity of the moving lattice. At low 
velocities, the sample is unaffected. Above some critical 
velocity Vc^ dissipation sets in abruptly. We determine 
Vc from a fit of Nc to the intersection of two lines with 
slopes and a: 



nd{v) 



nd{0) X [1 - max{0, a x {v — vc))] (1) 



The critical velocity which we obtain from this proce- 
dure is consistent for a large range of hold times, varying 
by less than 15% when the hold time t is changed by a 
factor of 20. We explore the BEC-BCS crossover by adi- 
abatically ramping the magnetic field to different values 
after evaporation and repeating the measurement as be- 
fore. The crossover is parameterized by the interaction 
parameter l//ci?a, where kp is the Fermi wavevector and 
a is the B-field dependent s-wave scattering length p^ . 
Again, we observe a threshold for dissipation. 

Figure [2] shows the measured critical velocity through- 
out the BEC-BCS crossover. The maximum near reso- 
nance is consistent with the picture of a crossover be- 
tween two different types of excitation, as discussed in 
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FIG. 2: Critical velocities throughout the BEC-BCS 
crossover. A pronounced maximum was found at resonance. 
Data is shown for a Vb = 0.2 Ep deep lattice, held for t=500 
ms. The solid line is a guide to the eye. 



the introduction, and confirms that superfluidity is most 
robust on resonance. 

To illuminate the role of the inhomogeneous density 
distribution, we performed experiments in which the en- 
tire sample was perturbed by a uniform lattice. Lattice 
beams with 80 /im waists probed a more tightly confined 
sample of 2 x 10^ pairs, with spatial extent Rtf — 37 jam. 
The onset of dissipation seen in Figure [3] is still striking, 
but now loss is observed at much lower lattice velocities, 
in spite of a larger Fermi energy Ep = hx 12 A kHz. More- 
over, the onset of dissipation is slightly more gradual. 
When the magnetic field was varied across the Feshbach 
resonance, we again found a maximum of the critical ve- 
locity near resonance. The lowering of the critical veloc- 
ity due to the inhomogeneous density profile is expected, 
since at lower density, both the speed of sound and (on 
the BCS side) the pairing energy decrease. Although the 
critical velocity should approach zero in the low density 
wings of the cloud, we still observe a sudden onset of 
dissipation at a finite velocity, similar to studies in Ref. 
[Sl, llO, UjJ , where a laser beam pierced through the whole 
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condensate, but in contrast to studies reported in 

In the limit of vanishing perturbation, the critical ve- 
locity should be given by the Landau criterion. In Fig- 
ure [4] we address the effects of a finite lattice potential in 
the original lattice configuration, as depicted in Figure 
[H The critical velocity is shown to be a decreasing func- 
tion of Vb, saturating in the limit of low lattice depth 
(Vb < 0.03 Ep)' This behavior is consistent with nu- 
merical simulation 0, [13] • Measurements at the small- 
est lattice depths had large uncertainties, as the hold 
time required to observe a heating effect of the lattice 
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FIG. 3: Effects of density inhomogeneity on the critical veloc- 
ity. A configuration in which the lattice beams (80/im) were 
larger than the trapped sample (37/im) results in loss in the 
condensate number Nc at significantly lower velocity. Data is 
shown for a Vb = 0.15 Ep deep optical lattice held for 200 ms 
at a magnetic field of 822 G. 



approached the natural lifetime of our sample. For this 
reason we studied the field dependence (Figure [2]) at an 
intermediate lattice depth, where Vc was more well de- 
fined. 

For comparison with theory we reference the local 
Fermi velocity at the trap center = Vp{l + P)'-^^^ 
— 39 mm/s, where v'^. = ^^EfItti is the Fermi ve- 
locity of a non-interacting gas at the trap center, and 
^ = —0.58 is a universal parameter characterizing uni- 
tarity limited interactions [itI. ITsI. [i9| . For vanishing lat- 
tice depth, the observed critical velocity at unitarity ap- 
proaches Vc/vf — 0.25. If we use the local Fermi ve- 
locity Vf^w at the waist of the lattice, we obtain 
'^c/'^F.w = 0-34 . The difference between these values 
indicates the uncertainty due to residual density inho- 
mogeneity. The local speed of sound in a Fermi gas at 
unitarity is 



Cs = Vf,ioc (1 + Py^yVs = 0.37 Vf,io 



The critical velocity for pair breaking is 
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0.34 v.^. 



(2) 



(3) 



with A = 0.50 vl^ij2m [13, lli and /i = (1 vl^iJ2m. 
These two values should provide approximate upper 
bounds to the critical velocity at unitarity [2, J^] . It seems 
natural that the combination of both excitation mecha- 
nisms lowers the critical velocity further. Within these 
uncertainties, and those of the density, the theoretical 
predictions agree with the experimental results. 

Up until now, we have deferred a discussion of how the 
moving lattice couples to the excitations. In a pure sys- 
tem at zero temperature, one would expect the excitation 
spectrum to exhibit discrete resonances, where the per- 
turbation couples only to modes with the /c- vector of the 
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FIG. 4: Critical velocities at different lattice depths. The 
results show Vc to be a decreasing function of lattice depth Vq. 
In the limit of low Vb, Vc converges to a maximum value of 0.25 
Vf- Data was taken near resonance, at 822 G {l/kpa — 0.15) 
for hold times t = 250 ms, 500 ms, 1000 ms, 2000 ms (squares, 
diamonds, circles, triangles). 



lattice. On the other hand, at finite temperature, it is 
possible that the lattice drags along thermal atoms which 
are point-like perturbations and can create excitations at 
all /c- vectors. Our observation that the dissipation sets 
in at a certain threshold velocity and increases monoton- 
ically with velocity is consistent with the participation of 
the thermal component. 

We further elucidated the role of thermal excitations, 
by varying the temperature. Gradually reducing the trap 
depth from U^io limine during exposure to a lattice mov- 
ing above Vc^ will suppress the accumulation of a thermal 
component. The lifetime in this case exceeded that for 
a sample held at a fixed depth of either /7o or Umin- 
For Bose-Einstein condensates, theoretical papers em- 
phasized the role of the thermal component in the Lan- 

[13,1211. This 



dau damping process in a moving lattice 
was confirmed qualitatively in an experiment at Florence 
(l2| in which the lifetime of the sample was drastically 
improved by eliminating the thermal atoms. 

In our experiments, the clouds heated up during the 
exposure to the moving lattice. Figure [5] shows the in- 
crease in the number of thermal atoms and the loss in 
the total number of atoms due to evaporative cooling. In 
an idealized model, where density is fixed, constant dis- 
sipation would result in a linear decrease in the number 
of atoms due to evaporative cooling. Our data show an 
accelerated decrease, possibly reflecting increased dissi- 
pation due to the increasing fraction of thermal atoms. 
However, an accurate model should include the change 
in density (and therefore critical velocity) during the ex- 
posure time. Additional impurity atoms (e.g. sodium 
atoms) could cause dissipation even at zero temperature 
and would allow more controlled studies of the dissipa- 
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FIG. 5: Number of pairs which remained in the condensate 
Nc (filled circles) and thermal component Nth (open circles) 
after being held in a Vb = 0.35 Ef deep optical lattice moving 
with velocity = 6 mm/s for a variable hold time. The ther- 
mal component shows a linear increase (dashed line), whereas 
Ncond showed an accelerated loss, and is fit to a quadratic 
function (soUd line). 

tion mechanism. Unpaired atoms in clouds with popula- 
tion imbalance may not play this role because of phase 
separation effects [22]. 

Another possible dissipation mechanism in a lattice is 
the creation of two excitations through a dynamical or 
modulational instability. Such an instability [23] occurs 
already for weakly interacting particles moving through 
a lattice with momentum when they collide and scat- 
ter into states with momenta q ± Sq^ analogous to opti- 
cal parametric generation [24]. This process is energet- 
ically possible only above 0.5 qs^ where qB = h/2\L is 
the Bragg momentum which defines the edge of the Bril- 
louin zone. This corresponds to a velocity v = 11 mm/s 
for fermion pairs (and twice this value for single atoms). 
Since the highest critical velocities we observe are slightly 
below this threshold, and strongly decrease already for 
relatively small Vq^ it is very unlikely, that dynamical 
instabilities play a role in our experiments. Moreover, 
such instabilities should be strongly modified by Pauli 
blocking. For our ratio of local Fermi momentum to 
the Bragg momentum of 0.9, the first band is nearly full 
in the center of the cloud. For Bose-Einstein conden- 
sates, it has been recently predicted [25] and experimen- 
tally shown [14] that strong interactions can lower the 
threshold for the dynamical instability, close to the Mott- 
insulator transition. The range of 1-D lattice depths ex- 
plored here {Vq < 2 Er) is far from the ID Mott-insulator 
regime. We have observed the loss of coherence which 



typically accompanies the superfiuid to Mott insulation 
transition to occur only beyond Vq 25 Er. 

In conclusion, we have used a novel optical lattice ge- 
ometry to determine critical velocities in the BEC-BCS 
crossover without the complications of strong density 
inhomogeneity. This configuration could be applied to 
studies in atomic Bose gases which so far have been lim- 
ited by the inhomogeneous density [g", T^, im, [H. In 
addition, it would be interesting to study dynamical in- 
stabilities for fermions and the role of Pauli blocking. 
The authors would like to thank Aviv Keshet for exper- 
imental assistance. This research has been supported by 
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